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We have studied the neutral-ionic transition in organic mixed-stack compounds. A unified model 
has been derived which, in limiting cases, is equivalent to the models proposed earlier, the donor- 
acceptor model and the ionic Hubbard model. Detailed numerical calculations have been performed 
on this unified model with the help of the density-matrix renormalization-group (DMRG) procedure 
calculating excitation gaps, ionicity, lattice site entropy, two-site entropy, and the dimer order 
parameter on long chains and the unified phase diagram has been determined. 
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I. INTRODUCTION 

Two physically seemingly different models have been 
proposed in the literature to describe the neutral-ionic 
(N-I) transition first observed in organic mixed-stack 
charge-transfer saltspl^- Assuming that the coupling be- 
tween the stacks is weak, the system is modelled as a 
linear chain in which two kinds of molecules alternate. 
Assuming furthermore, that it is sufficient to consider 
a single orbital per site, Torrance and Hubbard^ have 
suggested that the different nature of the two molecules 
can be attributed to different values of the on-site energy, 
while the on-site Coulomb repulsion U can be taken to be 
identical. Taking into account the finite transfer integral 
between neighboring molecules we are led to the so-called 
ionic Hubbard model described by the Hamiltonian 



U <A 



i i 

where cj^(cj^) is the usual creation (annihilation) oper- 
ator of electrons at site i with spin a, Uia — c\^c^^ and 
rii = riia is the occupation number at site i. 

When the number of electrons is exactly equal to the 
number of sites, the competition between the on-site en- 
ergy difference A and the Coulomb energy U will deter- 
mine whether the system is a band insulator (BI) or a 
correlated Mott insulator (Ml). It is easily seen by look- 
ing at the energy levels in the atomic limit, shown in 
Fig. m that if hopping processes can be neglected, it is 
energetically favorable to have doubly occupied odd sites 
and empty even sites when J7 < A (the energy of the pair 
is = J7 — A < 0), while in the opposite case [U > A) 
in the lowest-energy configuration every site is occupied 
by one electron (the energy of such a pair is £' = 0), 
but their spin can be arbitrarily oriented. In the band- 
insulator state both the charge and spin gaps are finite, 
while in the correlated Mott insulator only the charge 
gap is finite. The transition between these phases takes 
place at U — A. 

As was first pointed out by Fabrizio et al.^ for finite 
values of the hopping integral t the transition between 
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FIG. 1: Energy levels of the ionic Hubbard model in the 
atomic limit for U < A and U > A. 



these two states occurs in two steps. The charge gap 
closes at a critical value J7ci, but it reopens immediately, 
while the spin gap vanishes at a different value Uc2 > 
Uci- The transition at U = Uci is of Ising-like while 
the one at Uc2 is a Berezinskii-Kosterlitz-Thouless (BKT) 
transition. A dimerized bond-order (BO) phase exists 
between the two critical values. 

Since then the model has been studied in detail by sev- 
eral groups using both analytic and numeri cal methods. 
The most recent works are in Refs. |i|ll7"8'9"TolOIIl 
where further references can be found. After some con- 
troversy, by now consensus seems to emerge, the numer- 
ical resultsS'ii'i^ seem to support the picture with two 
transitions. The model has also been extended to include 
intersite Coulomb interaction^"' ^■^ with rather similar re- 
sults. 

In realistic charge-transfer salts, however, the intra- 
molecular Coulomb energy is presumably the largest en- 
ergy and it is reasonable to assume that its unique role 
is to forbid doubly ionized molecules. If that is so, the 
transition between the neutral and ionic phases is driven 
by other couplings. In the model studied in detail by 
Avignon et aZ.,^"* Girlando and Painelli,^^ and Horovitz 
and one of the authors^^ it is assumed that donor and 
acceptor molecules alternate along the chain. In its neu- 
tral state (D°), the highest occupied molecular orbital 
(HOMO) of the donor is filled by two electrons of oppo- 
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site spin. The ionization energy to remove an electron 
from this orbital and to create thereby a singly ionized 
donor is /. Taking into account the Coulomb repul- 
sion [/d the energy of the doubly ionized (D^+) state is 
2/ + [/d . The donor molecules can thus be described by 
the Hamiltonian 

donor sites ^^'j 
donor sites 

The energies are measured with respect to the neutral 
configuration. If Uu ^ /, its role is to forbid doubly 
ionized sites, empty donor orbitals. 

On the other hand the lowest unoccupied molecular 
orbital (LUMO) of the acceptor is empty in its neutral 
state (A°). The energy is lowered to —A, where A is the 
affinity, when the level becomes singly occupied, singly 
ionized (A^), while the doubly ionized (A^^) acceptors 
have a high energy due to the Coulomb repulsion Ua- 
The acceptor molecules are described by the Hamiltonian 

Ha = -A ^ ni + UA ^ ni^n^i . (3) 

acceptor sites acceptor sites 

If Ua ^ A, its role is to forbid doubly ionized sites, dou- 
bly occupied acceptor orbitals. The atomic energy levels 
of the donor and acceptor molecules are shown schemat- 
ically in Fig. [3 
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FIG. 2: Atomic energy levels of donor and acceptor molecules. 

Since / and A are fixed internal parameters of the 
molecules, and it is assumed that I— A > 0, the molecules 
are neutral in the ground state of the stack unless a tran- 
sition to an ionized state is driven by the intermolecu- 
lar Coulomb attraction V between the oppositely ionized 
donors (D+) and acceptors (A~), 

nG = -V (2-n,)n,+i 

donor sites , , , 

(4) 

acceptor sites 

In this atomic limit this transition should occur aX V — 
I- A. 

However, in addition to these terms one has to take 
into account the charge transfer between the donor and 



acceptor molecules described by 

HcT = t E + cI+K^C..) • (5) 

all sites, a 

Since for finite hopping the molecules are always at least 
partially ionized, this may smear out the sharp transition 
between the neutral and ionic phases. 

In fact, exact diagonalization on relatively short 
chains^^ and valence-bond calculationsi^iiS indicated 
that the transition remains of first order until t < 
0.15(7 — A). For larger t values a second-order or BKT- 
like transition was observed. The charge, spin and 
charge-transfer gaps seemed to vanish at the same criti- 
cal intermolecular Coulomb coupling Vc, but the charge 
gap reopens again. 

Not only the interaction driving the neutral-ionic tran- 
sition is different in the two models but also the char- 
acter of the transition seems to be different. The aim 
of the present work is to perform a more careful study 
of the neutral-ionic transition. First it is shown that a 
unified model can be derived which is identical to the 
donor-acceptor model in the limit when the intramolecu- 
lar Coulomb repulsion forbids doubly ionized donors and 
acceptors, and at the same time it is also a good approx- 
imation to the ionic Hubbard model when t ^ A, but 
U — \ can be arbitrary. Detailed numerical calculations 
are performed on this unified model with the help of the 
density-matrix renormalization-group (DMRG)i2i proce- 
dure calculating excitation gaps, ionicity, single-site en- 
tropy, two-site entropy, and the dimer order parameter on 
long chains and the unified phase diagram is determined. 

The setup of the paper is as follows. In Sec. II we trans- 
form the Hamiltonian of the two models into an effective 
spin-1 model and show their relationship. In Sec. HI we 
discuss how the quantum phase transition between the 
neutral and ionic phases can be best described and our 
numerical procedure is presented. The results are given 
in Sec. IV, and the conclusions are drawn in Sec. V. 



II. UNIFIED SPIN-1 HAMILTONIAN 

As mentioned earlier, in the donor-acceptor model, 
the on-site Coulomb couplings are large compared to the 
other characteristic energies in the problem, and the limit 
f^D, Ua 00 can be taken whereby the doubly ionized 
molecular states (empty donor and doubly occupied 
acceptor A^~) are forbidden. Since only three states per 
site survive, this model can be mapped onto an effective 
5=1 spin model. The appropriate mapping between 
the allowed electronic states and spin states is shown in 
Fig. El 

As is seen, the transfer of an electron with spin t or J. 
from the neutral donor at site 2j — 1 to the empty, neutral 
acceptor at site 2j, the process ^2j-i^2j ~^ ^2j-i^2j 
corresponds in the spin language to an exchange pro- 
cess |0)2,_J0)2, ^ |-l)2,„i|l)2, or \0),,_,\0),, ^ 
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FIG. 3: Mapping of the allowed electronic states of donor and 
acceptor molecules to S = 1 spin states. 



11)2 I~l)2j- '^^^ opposite processes are also possible. 
No hopping could, however, take place between a neutral 
donor and an ionized acceptor or an ionized donor and a 
neutral acceptor or between ionized neighbors when the 
spins of the electrons are parallel. 

However, one has to take into account that due to the 
fermionic nature of the electrons, these hopping processes 
appear with different signs. Acting by the charge-transfer 
term on a neutral donor-acceptor pair 

HcT 4j„i|4j_ix|0)2i-l|0)2i 

=t{4j-iAji ' 4j-u4,t)Io)2j-iIo)2j (6) 

= t(|T>2,-lU>2,-U)2,-l|T>2,). 

Similarly when the charge-transfer term acts on an 
^2P%+i pair 



HcT |0)2i4j + i|4i+u|0)2j + l 

=*(4jt4+u - 4ji4j+iT)lo)2j|o) 

= t(|T>2,U>2, + l-U>2,|T>2, + l). 



2j + l 



(7) 



The charge-transfer processes can thus be described in 
the spin language by the Hamiltonian 



all sites 



(8) 



all sites 



where and S[' are the usual raising and lowering oper- 
ators corresponding in the fermionic picture to removing 
or adding an electron. measures the state of the ion, 
and the product SfSf^i makes sure that no charge trans- 
fer takes place between a neutral and an ionized site. 

Using the commutation relations of the spin operators, 
this part of the Hamiltonian can be written in the form 



i 



(9) 



It is easy to check that the same matrix elements are 
obtained if the other terms of the Hamiltonian are writ- 



ten in the spin language as 

donor sites 

acceptor sites 
all sites 

Due to charge conservation 

donor sites acceptor sites 

and therefore 

all sites 



(10) 



(11) 



(12) 



We mention that the effective spin-1 Hamiltonian used 
in Ref. 16] differs somewhat from the one given above. 
The charge-transfer term was written in the form 



all sites 



(13) 



all sites 



The two expressions can be related by a unitary trans- 
formation, by introducing phase factors in the mapping, 
which alternate with a four-site periodicity. The appro- 
priate mapping is given in Fig. 01 
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FIG. 4: Phase factors in the mapping of the donor-acceptor 
model to a spin-1 model. 

In their work on the extended ionic Hubbard model 
Aligia and BatistciiS have shown that a similar mapping 
to an effective spin-1 model can be used in this model 
as well. When the number of electrons is exactly equal 
to the number of sites, the occupancy of states on even 
and odd sites is strongly correlated. Due to charge con- 
servation an empty odd site, which itself has low energy, 
implies that an even site is doubly occupied, and the en- 
ergy of this pair E = 2A -I- [/ is larger than that of a 
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pair of singly occupied even and odd sites {E = 0), or 
when the even site is empty and the odd site is doubly 
occupied {E = U — A). Near the expected transition, 
[/ « A, doubly occupied even sites and empty odd sites 
can, therefore, be neglected. 

With this provisio the electronic states of the ionic 
Hubbard model can be mapped to the three S — I states. 
The convention used is shown in Fig. [S] 
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FIG. 5: Mapping of the electronic states of the ionic Hubbard 
model to 5 = 1 spin states for the even and odd sites. 

The requirement that all matrix elements should be 
identical in the two representations leads to the following 
forms for the three terms of the Hamiltonian of the ionic 
Hubbard model: 



Ht = 



Hu = 



all sites 
all sites 

U (l-(^r)'), (15) 



odd sites 



.even sites 



odd sites 



.(16) 



The hopping term has exactly the same form as the 
charge-transfer term in the donor-acceptor model. The 
resemblance of the atomic term of the ionic Hubbard 
model to that of the donor-acceptor model becomes man- 
ifest when charge conservation is taken into account. The 
sum of the on-site Coulomb term (Tiu) and the ionic term 
(T^a) can be written as 



Hu + 'Ha 



where N is the number of sites and the parameter 

e = U -A 



(17) 



(18) 



has been introduced. Apart from a constant term this is 
the same as the atomic part of the donor-acceptor model 
given in Eq. (|12|l . if the identification 



6 = -{I -A) 



(19) 



is done. Thus in this representation the Hamiltonian of 
the donor-acceptor model is just an extension of that of 
the ionic Hubbard model, including the nearest-neighbor 
Coulomb interaction. Even sites in the ionic Hubbard 
model correspond to acceptors and odd sites to donors. 



There is, however, an essential physical difference. 
While in the ionic Hubbard model the two transitions oc- 
cur at positive values of e, in the donor-acceptor model 
s is assumed to be negative, and the nearest-neighbor 
Coulomb coupling is needed to drive the transition. The 
question we want to address in the remaining part of 
this paper is how the transitions obtained in the ionic 
Hubbard model can be related in this unified model to 
the first or second order transitions of the donor-acceptor 
model. 



III. DETECTING AND LOCATING PHASE 
TRANSITIONS 

A customary numerical procedure to detect and locate 
quantum phase transitions is to calculate energy gaps or 
correlation functions. A drawback of this procedure is 
that even if DMRG allows to determine these quanti- 
ties on long chains, quite often it is very difficult to draw 
firm conclusions from finite-size scaling. This is well illus- 
trated by the longstanding controversy over the existence 
of two phase transitions in the ionic Hubbard model. 

Recently we have proposed^^ that the two-site entropy, 
which is easily accessible in DMRG, should be considered 
when quantum phase transitions arc studied. In this sec- 
tion we present briefly the quantities that were used to 
locate the phase transitions and our method to control 
the numerical accuracy is also summarized. 



A. Energy gaps 

The natural quantities to be calculated are the excita- 
tion gaps. There are several types of excited states which 
can be considered. 

1. AEc is the energy needed to add or remove an elec- 
tron. Since originally there is an equal number of T 
and I electrons in the system, in the spin language, 
this gap is the energy difference between the lowest 
lying states with total spin — ±1 and 0. De- 
noting the lowest energy in the spin sector by 
EoiS^), AE, = EoiS^ = 1) - Eo{S^ = 0). 

2. AEg is the energy gap of spin-flip excitations. In 
the spin language, it corresponds to the energy gap 
between the lowest lying states with total spin 2 
and 0, A£:, = Eq{S^ = 2) - Eo{S^ = 0). 

3. AEcT is the excitation energy needed to trans- 
fer charge from a donor to an acceptor without 
changing the total charge and spin. In the spin 
language, it is the energy difference between the 
two lowest lying states of the ~ spin sec- 
tor, A£:cT = Ei{S^ = 0) - EoiS^ = 0). When 
calculating this quantity one has to take into ac- 
count that for an open chain with an even number 
of lattice sites the two lowest lying singlet states 
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are separated by a finite energy gap even in the 
thermodynamic limit, if the system is dimerized. 
In order to get the proper charge-transfer gap the 
finite-chain calculations have to be done with peri- 
odic boundary condition. 

In order to facilitate comparison with the work by 
Manmana et al^ where the behavior of the gaps of 
the nontruncated ionic Hubbard model have been care- 
fully studied, we note that their Ai is twice our charge 
gap, Ai = 2AEc, the spin gap is definied in the same 
way, As = AEg, while the quantity we characterized 
as charge-transfer gap is called by them exciton gap, 

Ae = Ai?cT- 

B. lonicity and dimer order 

Alternatively one can look at an appropriately defined 
order parameter. One such quantity could be the average 
charge on the molecules, the ionicity. In the spin repre- 
sentation it is given by the expectation value of (S'f)^, 

g,^{^Gs\iS!)'\^Gs), (20) 

where I^'gs) is the ground state wave function. It turns 
out that this is a good quantity to look at when the tran- 
sition is of first order and the ionicity has a finite jump at 
the transition point. In a second-order transition, how- 
ever, the ionicity is continuous, although with an infinite 
derivative in the thermodynamic limit, so it is difficult 
to get reliable information from its behavior in finite sys- 
tems. 

As has been pointed out in Refs. another proper 

order parameter could be the dimer order parameter D. 
An indication of the existence of dimer order can be ob- 
tained by measuring the alternation in the bond energy 

A = {'H,-l,^) - (^M+l) (21) 

in the middle of a long enough open chain. The dimer 
order parameter can thus be defined as 

D= hm \Dr,/2\. (22) 

N—^oo 

C. Single-site and two-site entropies 

Knowing the reduced density matrix pi of site i, which 
is obtained from the wave function of the total system by 
tracing out all configurations of all other sites, the von 
Neumann entropy of this site can be determined from 
Si = —Tt Pi In Pi. For a model with q degrees of free- 
dom per site Si may vary between and Inq. In a chain 
with free ends, the site entropy at the center of the chain 
Sjv/2 is free from end effects if the chain length is larger 
than the coherence length, and the anomalies appearing 



in this quantity as a function of the coupling constant 
can be used to detect quantum phase transitions. The 
site entropy has a jump at a first-order transition, while 
if a second-order transition takes place between two dif- 
ferently ordered phases, the site entropy is expected to 
show a sharp maximum at the transition point. A simi- 
lar quantity, the concurrence, can also be used to locate 
quantum phase transition points, as has been shown by 
Vidal et al}^ for the Ising model in transverse field. 

The single-site entropy, however, is not always a good 
indicator of a quantum phase transition. It can be fea- 
tureless when it is insensitive to the breaking of symme- 
try that distinguishes the two phases. Recently it has 
been pointed outi^ that in such cases quantum phase 
transitions can be more efhciently detected and located 
by studying the von Neumann entropy of an ensemble 
of two lattice sites in the middle of a long chain. The 
quantity defined by 

Si,i+i = -Tipis+i \npis+i , (23) 

where is the reduced density matrix of two neigh- 

boring sites, may display — when taken in the middle of 
the chain, at i = N/2 — a relatively sharp maximum at a 
phase transition point even though the single-site entropy 
does not. 

When the system gets dimerized at the transition, 
the breakdown of translational symmetry can also be 
detected by calculating — as an alternative to the usual 
dimer order parameter — the difference of two-site en- 
tropies on neighboring sites, in the center of the chain 

Ds = s»+i,,;+2 - i^N/2. (24) 

D. Numerical accuracy 

The numerical calculations were performed on finite 
spin chains up to iV = 800 lattice sites when open bound- 
ary condition (OBC) was used and up to = 256 
sites for periodic boundary conditions (PBC) using the 
DMRG techniqucr^ and the dynamic block-state selec- 
tion (DBSS) approachi22i2i All data shown in the figures 
were obtained with OBC unless stated otherwise. We 
have set the threshold value of the quantum information 
loss X to 10"^ and the minimum number of block states 
Afmin to 256. The maximum number of block states 
varied in the range 300 - 800 for OBC and 500 - 2000 
for PBC, respectively. All eigenstates of the model have 
been targeted independently using two or three DMRG 
sweeps. The entropy sum rule was checked for all finite 
chain lengths for each DMRG sweep, and it was found 
that the sum rule was satisfied after the second sweep 
already. 

After accomplishing the infinite lattice procedure and 
using White's wave-function transformation method— 
the largest value of the fidelity error of the starting vec- 
tor of the superblock diagonalization procedure, Se^s,^^^ = 
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1 — (^'xl^'stv), where ^stv is the starting vector and 
is the target state determined by the diagonahzation of 
the superblock Hamiltonian, was of the order of 10^^. 

As another test of the accuracy we calculated the ionic- 
ity (gi), lattice site entropy (si), two-site entropy (si^i+i) 
and dimer order (-Di) profile using PBC. In principle, 
each of them should be independent of i for a given finite 
chain length. Using the DBSS approach with x = 10^^, 
Mjnin = 256 for chains up to iV = 256 sites the values 
obtained differ from their mean by less than 10~^ for all 
sites. 

It is important to emphasize that when a maximum 
or cusp is looked for in the single-site or two-site en- 
tropies, the application of the DBSS approach is crucial, 
otherwise — if the number of block states m were fixed 
during the calculation — a very large number of states 
would have to be used. This is due to the rapid increase 
of entanglement around the transition points, thus a con- 
stant m can lead to a significant and non-constant cut in 
the entropy functions. 

In order to obtain the energy gaps, order parameters 
and entropies in the N ^ oo thermodynamic limit, finite- 
size scaling analysis has been performed using scaling 
functions appropriate for OBC and PBCi^ 



IV. NUMERICAL RESULTS 

Since a nontrivial truncation procedure has to be per- 
formed to reduce the ionic Hubbard model to the unified 
model with e > and V = 0, first we have studied to 
what extent the truncation modifies the location of the 
phase transition in the model. As a next step we checked 
the first-order transition of the donor-acceptor model for 
e < and t V. Then we extended the calculation to 
intermediate values of the couplings to see how the two 
limits can be related. 



A. The ionic Hubbard model versus the unified 
model at = 

As mentioned before the ionic Hubbard model given 
by Eq. ^ has been shown to possess two phase transi- 
tions. An ionic bond-ordered dimerized phase has been 
found between the neutral regular and the ionic regular 
phases. Since the location of the BKT transition is an 
especially difficult problem when the vanishing of gaps 
is studied numerically, we have chosen to calculate the 
single-site and two-site entropies as a function of U for 
various values of t by taking A as the unit of energy. As 
an example we show in Fig.^lour results for various finite 
chain lengths at t/A = 0.5. 

As seen in the figure the single-site entropy of the cen- 
tral site is a rather smooth curve without any cusp or 
sharp maximum, while the two-site entropy has two max- 
ima. As the size of the system is increased the maxima 
get closer, but finite-size scaling analysis shows that they 
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FIG. 6: Single-site and two-site entropies as well as the dimer- 
ization in the two-site entropy of the ionic Hubbard model for 
t/A = 0.5 as a function of U/A. 



remain separated and develop into two cusps in the ther- 
modynamic limit. The same behavior is found for other 
values of t/A, except that with decreasing t the two peaks 
get even closer. This result is in perfect agreement with 
the two-transition scenario mentioned above. One peak 
can be identified with Ud, the other with Uc2- The phase 
diagram of the ionic Hubbard model obtained from these 
calculations for t/A < 0.5 is given in Fig. [T] 

The values obtained by Manmana et al.^ for Uci and 
the bounds for Uc2 from their study of the gaps are also 
shown in the figure. They fit very nicely on the curves 
obtained from the peaks of the two-site entropy. 

A hint about the nature of the phases can be obtained 
when the alternation of the two-site entropy is looked 
at. As seen in the third panel of Fig. Elit vanishes for 
small and large values of U /A and is finite in the in- 
termediate region only, indicating that in the BO phase, 
the system is spontaneously dimerized, while it is regu- 
lar in the BI and MI phases. The critical value where 
dimerization first appears coincides well with Ud, while 
it cannot be established with certainty that the dimer 
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FIG. 7: Phase diagram of the ionic Hubbard modeL The 
critical values at t/ A = 0.05 and 0.25 are taken from Ref. 0. 



order disappears exactly at Uc2, since convergence at the 
BKT transition is very slow. There is no doubt, how- 
ever, that the dimer order disappears at large U values 
and a consistent picture is obtained if it is assumed that 
spontaneous dimerization occurs in a narrow range only, 
between the two transitions. 

Furthermore, the phase boundaries allow to determine 
the range where the truncated Hamiltonian can be jus- 
tified. That is the regime where (C//A)ci and ([//A)c2 
vary linearly with t/A, i.e., where the relevant parameter 
is the combination e/t = {U — A)/t. The slopes of the 
curves give 



1.4. 



(25) 



and the curvature is rather small up to t/A w 0.5. 

As a check of these results we have repeated the cal- 
culations on the unified model at = for e > 0. The 
results for various finite chain lengths are shown in Fig.|Sl 

The site entropy is a rather smooth curve, like for 
the full model. Even if eventually it develops a cusp 
at {e/t)ci — 1-3, there is no indication in the single-site 
entropy of the second, BKT-like transition. In contrast 
to this, the two-site entropy exhibits again two maxima. 
Although for the longest chains the entropy curves have 
been smoothened close to the maxima due to the fact 
that the number of block states selected dynamically has 
reached the maximum number of block states that our 
code could handle, nevertheless finite-size scaling anal- 
ysis of the curves gave {e/t)ci « 1.3 and (e/i)c2 ~ 2 
for the two critical values, in agreement with our result 
given in Eq. (|25|l and with earlier results on the full ionic 
Hubbard model>ii*i^ 

Moreover, the two-site entropy Ds shows clearly that 
the system is spontaneously dimerized above the first 
critical value (e/t)ci, but again although Dg certainly 
vanishes at large enough e/t, the exact location where 
this happens cannot be determined from the available 
chain lengths using finite-size scaling analysis. 
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FIG. 8: Same as Fig. |S|but for the effective spin-1 model at 
= as a function of e/t = ([/ - A) /t. 



B. Large negative values of e/t 

The donor-acceptor model corresponds to e < 0. In 
the absence of charge transfer between neighbors [t = 0) 
the transition from the neutral to the ionic phase takes 
place at 



V=Ul-A) 



(26) 



This is a first-order transition where the ionicity jumps 
from zero to one and none of the gaps go to zero con- 
tinuously as discussed in Ref. In the neutral state, 
all sites are in a pure state, the single-site and two-site 
entropies are zero. In the ionic phase the high degener- 
acy is due to the arbitrary orientation of the spins, and 
the site entropy is In 2, the two-site entropy is In 4. The 
jump occurs at the transition point. 

For finite t the ionicity and the single-site and two-site 
entropies are nonvanishing in the neutral phase and their 
value is less than the maximal one in the ionic phase. 
Nevertheless, when t is small enough compared to V and 
e|, a finite jump can be detected in these quantities near 
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V w — e/2. The behavior of these quantities as a function 
of V/t for e/t — —10 is shown in Fig.|51 
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FIG. 9: The single-site entropy (si), two-site entropy (si,i+i), 
and ionicity {gi) in the middle of the chains for different chain 
lengths at e/f = — 10 as a function of V/t. 

The excitation energies obtained on finite chains using 
PBC for e/i = — 10 as a function of V/t in the vicinity 
of the N-I transition are shown in Fig. ^| Although 
the charge gap AE^ has a minimum for all finite TV, it 
remains finite even in the thermodynamic limit. The 
minimum in the N ^ oo limit is at V/t ~ 5.13, the 
same value where the ionicity and the entropies have a 
jump. At the same point, the spin gap AEs and the 
charge-transfer gap Ai^cT jump to zero and remain zero 
for larger V/t values. 

The dimer order parameter was found to be of the 
order of 10^^ or less for aU V/t for N = 400 lattice sites 
already. Finite-size scaling extrapolation gives less than 
lO^'* for the dimer order parameter. 

Similar results have been found everywhere for e/t < 
—8 except that the transition point shifts to larger V/t 
values with decreasing e/t as expected. We can, there- 
fore, conclude — in agreement with earlier worksi^*i& — 
that for small values of t/V, or equivalently for large 






FIG. 10: Charge gap, spin gap and charge-transfer gap for a 
small nonvanishing charge-transfer integral {e/t = —10) using 
PBC. 



negative values of e/t the system undergoes a first-order 
transition from a regular neutral to a regular ionic phase. 



C. Smaller values of e/t 

As we have seen, the unified model describes on the 
one hand the two transitions taking place at V = and 
e/t > 0, and on the other hand, the unique first-order 
transition at large negative values of e/t. The question 
we want to address now is how these two limits can be in- 
corporated into a unified phase diagram. We have, there- 
fore, done similar calculations for fixed smaller negative 
and positive values of e/i as a function of V/t. The exci- 
tation energies obtained on finite chains using OBC for 
e/t — —2 are shown in Fig. 1111 

Comparison to the results obtained in Ref. Q shows 
that the gaps behave surprisingly similarly to that found 
in the full ionic Hubbard model. The charge-transfer gap 
vanishes at a critical value {V/t)ci which for e/t — —2 
is {V/t)ci ~ 1.39, but reopens again and closes a second 
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FIG. 11: Charge gap, spin gap and charge-transfer gap for 
e/t = -2. 



has to be studied, as has been done in the ionic Hub- 
bard model. The behavior is again very similar. The 
extrapolated value of the measured dimer order parame- 
ter is exceedingly small, of the order of 10~^ — 10~^ for 
V/t < {V/t)ci. It becomes finite above this transition 
point, but disappears again at a larger value of V/t. Cal- 
culations on very long chains {N = 800) still do not allow 
a reliable extrapolation to determine the coupling where 
this happens. 

Since it is difficult to locate the second transition from 
the vanishing of the spin gap or the dimcrization, we have 
studied the behavior of the entropy. Our results obtained 
on finite chains are shown in Fig. 1121 



^ 40 
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time at another critical point, whose exact location, how- 
ever, cannot be determined by finite-size scaling from the 
available chain lengths. When the same calculations are 
performed with PBC, the charge-transfer gap vanishes at 
{V/t)ci and remains zero. Thus the reopening of the gap 
is due to our use of OBC. As we will see, the system is 
spontaneously dimerized for V/t larger than this critical 
value, and as has been discussed in Sec. Ill, this gives 
rise to the finite gap between the two lowest lying levels. 

The extrapolated spin gap, on the other hand, is finite 
at {V/t)ci and closes at a second critical value {V/t)c2- 
Due to the very slow decay of the spin gap even on rather 
long chains, the exact location of the second phase tran- 
sition could not be determined accurately from the gap. 

The charge gap AEc has a minimum for all N and the 
location of the minimum scales in the thermodynamic 
limit to {V/t)ci — 1.39, however, the gap remains finite 
at this transition point. The same behavior has been 
found for the one-particle gap Ai in Ref. 1^. 

The ionicity (not shown) is no longer a good indicator 
of the neutral-ionic transition, it is found to be a con- 
tinuous function of V/t. Instead of that the dimer order 




FIG. 12; The single-site entropy, two-site entropy, and the 
dimerization of the two-site entropy for e/t = —2. 

The single-site entropy seems to develop a cusp indi- 
cating a second-order phase transition taking place at 
{V/t)ci — 1.4, but only one transition is found. In con- 
trast to this the two-site entropy possesses two maxima 
for all finite chain lengths and they both develop into a 
cusp-like peak in the thermodynamic limit. The criti- 
cal points are at V/t ~ 1.39 and 1.66. Calculations in 
the vicinity of V/t = 1.39 on chains up to = 800 lat- 
tice sites confirm that although the rise in the two-site 



10 



entropy is rather steep, there is no jump in it, the tran- 
sition is of second order. 

Similar results, two transitions have been found for 
other values of e/t as a function of V/t, when e/t is 
larger than about —8. When e/t is larger than the crit- 
ical {e/t)ci « 1.3 or {e/t)c2 ~ 2, one or both transitions 
appear for negative V/t values. 

D. Phase diagram 

The results obtained until now can be summarized in 
the phase diagram shown in Fig. 1131 
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Second order transition 
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Tricritical point 
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FIG. 13: Phase diagram of the unified model describing the 
ionic Hubbard and donor-acceptor models as a function of 
the parameters e/t and V/t. Dashed line indicates first-order 
transition, solid lines correspond to second-order or BKT 
transitions. The star denotes the tricritical point. The phase 
boundaries of the full ionic Hubbard model are also shown. 

The two transitions obtained in the ionic Hubbard 
model and the first-order transition of the donor-acceptor 
model appear in fact as two limiting cases of the unified 
model. When e is negative and |e| ^ t, the neutral-ionic 
transition is discontinuous, the ionicity jumps by a finite 
amount. A tricritical point appears at about —e/t « 8 
and V/t w 4.2, beyond which two transitions are found. 
The charge-transfer gap vanishes at one of the transi- 
tions, while the spin gap does so at the other one. 

Since the dimer order parameter is nonvanishing in the 



narrow region only, between the two continuous transi- 
tions, three phases can be identified, a regular neutral, a 
regular ionic, and a dimerized phase inbetween. 



V. CONCLUSION 

In summary, we have studied the neutral-ionic transi- 
tion in organic mixed-stack compounds using the density- 
matrix renormalization-group method. First, we have 
shown that a unified model can be derived which is iden- 
tical to the donor-acceptor model when the intramolecu- 
lar Coulomb repulsion is the largest energy in the prob- 
lem and it represents a good approximation to the ionic 
Hubbard model when A/t ^ 1. 

Detailed numerical calculations have been performed 
on this unified model calculating excitation gaps, ionic- 
ity, lattice site entropy and dimer order parameter on 
long chains. We have shown that the best quantity to 
study is the two-site entropy. It exhibits the transitions 
most clearly, even at points where finite-size effects are 
important and the vanishing of gaps in the thermody- 
namic limit is difficult to establish. In this way we could 
determine the unified phase diagram as a function of in- 
tersite Coulomb interaction and the intraatomic energies. 
The results are in complete agreement with earlier results 
on the ionic Hubbard model, while earlier works on the 
donor-acceptor model have found a single (first or second 
order) transition only. This is probably due to the close- 
ness of the two transitions, which could not be resolved 
on short chains. 

In the model studied, spontaneous dimerization ap- 
pears in a narrow range of couplings only, while exper- 
iments indicate that the ionic phase is dimerized when 
electron-phonon coupling is taken into account and the 
displacement of ions is permitted. Extension of the cal- 
culations in this direction is in progress. 
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